Diffusion in the general theory of relativity 
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The Markovian diffusion theory in the phase space is generalized within the framework of the gen- 
eral theory of relativity. The introduction of moving orthonormal frame vectors both for the position 
as well the velocity space enables to bypass difhculties in the general relativistic stochastic calculus. 
The general relativistic Kramers equation in the phase space is derived both in the parametrization 
of phase space proper time and the coordinate time. The transformation of the obtained diffusion 
equation under hypersurface-preserving coordinate transformations is analyzed and diffusion in the 
expanding universe is studied. It is shown that the validity of the fluctuation-dissipation theorem 
ensures that in the quasi-steady state regime the result of the derived diffusion equation is consistent 
with the kinetic theory in thermodynamic equilibrium. 
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1. INTRODUCTION 

The theory of Markovian diffusion within the frame- 
work of the relativity theory has re-attracted consider- 
able interest during the last decade. Compared with the 
non-relativistic case, relativistic Markovian diffusion be- 
comes significantly more difficulty due to conceptual and 
technical problems. Over the years many pap ers have 
been devoted to this issue (see e.g. [l| - [12]), a few 
studied also diffusion in gravitational fields within the 
theory of general relativity - However, at least 
up to the knowledge of the author, a general accepted 
formulation of the theory of Markovian diffusion within 
the framework of general relativity is still missing. Most 
investigations of statistical processes in relativistic fluids 
in cosmology and astrophysics used alternative physical 
approaches given by relativistic kinetic theory or phe- 
nomenological relativistic thermodynamics 13|-[l^. In 



particular the role of dissipative processes in the early 
evolution of the Universe has been studied by using non- 
equilibrium thermodynamics [23|-|2l|- The relativistic 
Boltzmann equation is an integral-differential equation 
which is difficulty to solve in the non-equilibrium case far 
from equilibrium. On the other hand, in the probabilis- 
tic approach Fokker-P lack- type equations are differential 
equations which can be much simpler solved. Besides, 
the probabilistic theory of diffusion processes within the 
framework of general relativity exhibit a fundamental in- 
terest and play a significant role in astrophysical and cos- 
mological problems (see e.g. ji^l - ). 

A crucial factor in relativistic diffusion is the fact that 
the velocity space in relativity is a hyperboloid (or a 
special 3-dimensional Riemannian manifold) embedded 
into the 4-dimensional velocity Minkowski space. In the 
derivation of the relativistic diffusion equation this spe- 
cific feature has to be taken into account in appropriate 
way, otherwise inconsistent results arise. One of the main 
difficulties in the derivation of a consistent Markovian dif- 
fusion equation in relativity is to handle the fundamen- 
tal Wiener process in a stochastic differential equation 



on non-Euclidian manifolds in a rigorous way. There ex- 
ist a well developed and rigorous mathematical theory of 
stochastic differential equations and diffusion processes 
in the base space on Riemannian manifolds with a defi- 
nite metric signature [l^, [s^]- However, this stochastic 
calculus can not be applied for manifolds with indefinite 
metric. Recently, the author presented a physically mo- 
tivated modification of this calculus to describe diffusion 
within the framework of special relativity in the the phase 
space of position and velocity. [3l| . In the present paper 
this formalism will be extended to a theory within the 
framework of general relativity in external gravitational 
fields. The main aim of the paper is the derivation of 
a generalized Kramers equation in external gravitational 
fields within the framework of general relativity theory 
both in the parametrization of the phase-space proper 
time and of the observer time within the approach of 
spacetime decomposition. The transformation property 
with respect of foliation-preserving coordinate transfor- 
mations is analyzed and diffusion in the evolving universe 
is studied. It is shown that the quasi-steady solution is 
in agreement with the kinetic theory in thermodynamic 
equilibrium. 



The paper is organized as follows. In Sec. 2 the princi- 
pal results pertaining to the special relativistic diffusion 
process is briefiy described. In Sec. 3 the formalism 
of Markovian diffusion in gravitational fields within the 
framework of general relativity is presented both in the 
parametrization of the phase-space proper time and the 
observer time. In Sec. 4 the transformation property of 
the derived diffusion equation is studied. As an example 
for the application of the formalism diffusion processes 
in the expanding universe are discussed in Sec. 5. and 
the conclusions are presented in Sec. 6. 
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2. THE SPECIAL RELATIVISTIC DIFFUSION 
PROCESS 

Let us first review the principal results pertaining to 
the special relativistic case 3l| used as the basic for the 
generalization within the frame of general relativity. 

The important point in relativistic diffusion is the 
observation that the velocity space in special relativ- 
ity is a noncompact hyperbolic 3-dimensional Rieman- 
nian manifold embedded into the 4-dimensional velocity 
Minkowski space. Using normalized velocity variables 
{fi — 0, 1, 2, 3) the hyperbolic metric structure for a rel- 
ativistic system of massive particles is described by the 
relation 



(^0)2 _ (^1)2 _ (^2)2 _ (^3)2 



1. 



(1) 



Therefore relativistic Markovian diffusion processes can 
be described in a rigorous way by using the mathemat- 
ical stochastic calculus on Riemannian manifolds, but 
adopted to the velocity space. 

Stochastic differential equations in diffusion theory 
with continuous pathway are defined by the fundamen- 
tal Wiener process W''(t). On a Riemannian manifold 
the fundamental Wiener process is difficult to handle. 
The key idea in the mathematical concept of diffusion on 
general Riemannian manifolds with definite metric signa- 
ture is to define a stochastic process on the curved man- 
ifold using the fundamental Wiener process each com- 
ponent of which is a process in the Euclidian space R'' 
[29j,(30j. The tangent space of a Riemannian manifold 
is endowed with an Euclidian structure and therefore we 
can move the manifold in the tangent space by construc- 
tion of a parallel translation along the stochastic curve 
with the help of the orthonormal frame vectors = 
ejj(x)9i (i,a=l...d) and the Christoffel connection coeffi- 
cients r^jj, X = {xi...Xd),di = d/dx^. In local coordinates 
on a Riemannian manifold the infinitesimal motion of a 
smooth curve c* [t) in M'' is that of 7* {t) in the tangent 
space by using a parallel transformation: dd = e^(x)d7'' 
and (ie^(x) — — T^je^jdc™. Therefore a random curve 
in the stochastic mathematical calculus on Riemannian 
manifolds in the position space can be defined in the same 
way by using the canonical realization of a d-dimensional 
Wiener process (defined in the Euclidian space) and sub- 
stituting ^7" — > dW°'{t). Thus the stochastic differential 
equations describing diffusion on a Riemannian manifold 
in the orthonormal frame bundle 0{M) with coordinates 
0{M) ={x%ej,} are given by [ii,^ 

dx\T) = e'^{T)odW^ +h\T)dT, (2) 

deL(r) = -r^ie\odx"^{T). 

Here r^j,(x(T))e^,(x(T)) = g^^ ,d,ei = -Y\^e\, g^^ is 
the Riemannian metric and 5°'^ the flat Euclidian metric 
where S""^ is the Kronecker symbol. 



The general mathematical model given by Equ.(2) de- 
scribes diffusion in the position space. It can not be 
generalized to pseudo-Riemannian manifolds in the gen- 
eral relativity theory. But Equ.(2) exhibit a significant 
difference to the non-relativistic Langevin equation in 
physics, where the stochastic force acts directly only on 
the change of the velocity and not on the position co- 
ordinates. Therefore the mathematical stochastic cal- 
culus for Riemannian manifolds can be applied with 
an appropriate modification adopting it to the velocity 
space. This requires the introduction of a moving ve- 
locity frame EI{t). A relativistic generalization of the 
Langevin equations can be defined in the fiber bundle 
space F{Ml) = {x\u\ E^} hy ^ 



dx^ir) — u^T)dT, 

du'ir) = El{T)odW''{T)+F\T)dT, 

dEl{r) - -7™i(u)i?i(r)od«"(r). 



(3) 



Here r is an evolution parameter along the world lines 
of the particles which can be chosen as the proper time. 
The laboratory time t = tu^ /c is a function of the 
proper time r and vP which here and below is defined 
by u° = [I + {u^Y + [u^f + [u^Y]^. 7j„,(u) are the 
Christoffel connection coefficients on the hyperboloid and 
F^=K^ /m, where are the spatial components of the 
4-force {i — 1,2,3), m is the rest mass of the particles 
and the a, b numbers the moving frame vectors Ea in the 
hyperbolic velocity space {a,b — 1,2,3). The stochas- 
tic force is described by the fundamental Wiener process 
with {W) = and the correlator {W''{t)W''{t + s)) = 
DsSab, is defined by an empirical diffusion constant D 
which is independent on the velocity. Stochastic integrals 
related with Equ.(3) are defined in the Stratonovich cal- 
culus denoted by the symbol o. Since the manifold on 
the hyperboloid is embedded into the Minkowski space, 
metric Gij{u) and connection coefficients 7j7;(u) in the 
velocity space are given by 



G.,(u) 
7j.(u) 



(4) 



Associated with the diffusion process described by 
Equ.(3) is a diffusion generator A.p(^]^,jj^^ 



D 



(5) 



where the fundamental horizontal vector fields Ha and 
on the fiber bundle F{M]^) are given by 



(6) 



We project the stochastic curve from the fiber space 
F{Ml) with coordinates r ={a;^u%£'^} to the phase 
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space with coordinates {x^,u^} : AF{M^)f{r) = 
Ap/(x, u, 0), where the diffusion generator in the phase 
space Ap is given by 



A, 



The generator Ap describes how the expected value 
(P{t,x) = E^lfCKr)] of any smooth function / evolve in 
time and satisfies the following equation: 



d_ 



(P(t,x,u) = Av3(r,x,u), 



(8) 



with (p(0,u)=/(x). Equ. (8) is a Kolmogorov backward 
equation. A Fokker-Planck equation (or forward Kol- 
mogorov equation) describes how the probability density 
function 0(r,x) evolves with time and is determined by 
the adjoint of the diffusion operator Ap. Therefore the 
special relativistic diffusion equation in the phase space 
takes the form 3l|: 



90 



div.AFa 



D 



-A,, 



(9) 



where A„ is the Laplace Beltrami Operator of the hyper- 
bolic velocity space given by 



depending not only on the velocities v'^ but also on the co- 
ordinates x'' . The complication arising by this fact can be 
bypassed by using the orthonormal frame vectors cm = 
e^^{x)d^ in the position space, where the subscript M, N 
numbers the vectors {M,N = 0,1,2,3), /i their compo- 
nents in the coordinate basis = 5^ = 0, 1, 2, 3) and 

tion 



The frame vectors satisfy the condi- 



,?^«^V^t,(x)e5^(x) 
9nv^^{x)e.)^{x) 



(13) 



where ?7jvfjv — diag{— 1,1,1,1) is the metric of the 
Minkowski space. The dual basis of the frame fields 
e^v/ are cotangent frame 1-forms 9 
the orthogonality relation 



M 



off dx^^ satisfying 



l,{x)9t'(x)=6^^. 



Each vector refer to the coordinate basis e,. 



(14) 



dfj, is as- 



signed a vector refer to the frame basis cm according to 
the rule 



{x)v 



M 



(15) 



Correspondingly, we find for v 



M 



du^dui 



^ 1 ^(VGG'^^), 



(10) 



and the corresponding divergence operator is given by 



(16) 



By using the orthogonal moving frames of the pseudo- 
Riemannian manifold in the position space the relation 
(12) gets the same form as in the special- relativistic case 
in Equ. (1): 



with G = det{Gy } = {uq)-"^ and G'^ = S'^ + v^v? . 



g^,{x)v^^v- = g^,,{xyii{xY^{x)v'^v'' (17) 
(11) = ^mnv'-'v'' = 1. 

The covariant derivative in the natural frame is defined 
as usual. 



3. GENERAL RELATIVISTIC DIFFUSION 
EQUATION 

The special relativistic formalism briefiy described in 
section 2 can be extended to a theory within the frame- 
work of general relativity. Besides orthonormal frame 
vectors in the velocity space in the generalization of his 
formalism the introduction of orthonormal frame vectors 
in the position space is required. 

In the case of general relativity the four- velocity of 
massive particles in the presence of a gravitational field 
with a metric gy,v{,x) satisfies the condition 

g^,,{x)v^v'' = 1, (12) 



VX=5,.i;''-fr>'', (18) 

where are the Christoffel connection coefficients (or 
coordinate connections). The covariant derivative of vec- 
tors refer to the frame basis can be written as: 

vX' = aX' + f]>^, (19) 

where fi*^ are called coefficients of spin connection (or 
frame connection). The relation between these two kind 
of connections arise from the metric compatibility condi- 
tion, which here can be expressed by V^O^ = 0, V^eJJ^ = 
0. From this condition one gets [321, 1S3| 

^f.Ni^) = C(x)r;:^e^(x) + el\x)dMx). (20) 
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The parallel transport of the components of the 4- 



vector V in the moving frame ba and the line elements 
dxf^ are given by 



(21) 



3';^j{x)v^Ut. 



As seen in Equ.(21) the change of the spatial compo- 
nents of the velocity of a particle is determined by 
the gravitational force described by the spin connection 
coefficients rt™j^{x). In addition one has to take into ac- 
count that the velocity is changed by a stochastic force 
^noise driven by the Wiener process In order to 

avoid difficulties in the description of the Wiener pro- 
cess on the velocity hyperbolic manifold the mathemat- 
ical stochastic calculus on Riemannian manifolds has to 
be applied where the Wiener process is moved along the 
orthonormal frames E^{u) [3l| in the three-dimensional 
hyperbolic velocity space defined by the relation 



or equivalently 



o=l 



Sab 



(22) 



(23) 



where Gij is the Riemannian metric of the hyperbolic ve- 
locity space and G*^ the inverse matrix of Gij . The met- 
ric and the Christoffel coefficients in the velocity space 
are defined in Equ. (4). The infinitesimal motion of the 
velocity v™{t) in phase space can be described by that 
of u°'{t) (a — 1, 2, 3) in the moving velocity frame E"^ by 
using the parallel transformation 



dE^^ir) = -Y:MEl{T)dv^' 



(24) 



A random curve in the phase-space can be defined as 
in section 2 by using the Wiener process substituting 
— >■ dW'^{T). Therefore analogous as in Equ. (3) in 
a consistent description of Markovian diffusion in gen- 
eral relativity the noise force is described by F™^^^ = 
EJ^{t) o dW"". A remarkable feature of Markovian diffu- 
sion on a Riemannian manifold is the supposition that 
for the diffusion coefficients only the orthonormal frame 
coefficients E^ are admissible which are directly related 
to the hyperbolic geometry of the velocity space. In con- 
trast, on Euclidian manifolds a much more general class 
of diffusion coefficients are permitted. Therefore gener- 
alizing Equ. (3) the stochastic differential equations de- 
scribing diffusion in a gravitational and external force 
fields (general relativistic Langevin equations) are given 
by: 



dx' = e\i{x)v^^dT, (25) 
dv" = E'^{T)odW -n';^j^{x)elj{x)v^v^'dT + 

dEZ\T) = --^Z{^,)Ei{r)odv-. 

Here an additional external force FI^^=Kl^^/m is taken 
into account where K^,^. are the components of the ex- 
ternal 4-force in the coordinate frame and m is the rest 
mass of the particles. In the moving frame this force is 
expressed by = O^Fj^^. The Christoffel connection 
coefficients on the hyperboloid Yni i^) defined in Equ. 
(4) and r is a parameter defined along the world line of 
the particles, which can be chosen as the phase-space 
proper time. 

Sufficient and necessary conditions for the existence 
and uniqueness of the stochastic differential equation (25) 
are that the drift and diffusion coefficients satisfy the 
uniform Lipschitz condition and the stochastic process 
X(T)={x(r), v(t)} is adapted to the Wiener process 
W'^(r) , that is, the output X(t2) is a function of W°(ri) 
up to that time (ti < T2). 



The diffusion operator A 



F{M) 



in the fiber bun- 



dle F{Ml) with coordinates r u™, for the 

stochastic process described by Equ. (25) is defined as in 
Equ. (5) with horizontal vector fields Ha and Hq derived 
analogous as in Sec. 2: 



= e: 



d 
dv" 



d 



dEV 



(26) 



d 



d 



lnli^)K F^- 



dE- 



The operator A^ja/) can be projected to the phase 
space with coordinates r —{x^,v™-} by A^(Af)/(r) = 
Ap/(x, u), where the diffusion generator in the phase 
space Ap is given by 



Ap - ^^'i;: 



d 



d 



d 



dv 



dx^ 



(27) 



Here the first term contains the Laplace-Beltrami op- 
erator A„ = S''''E^g^EJ^-^ in the hyperbolic ve- 
locity space. The force is composed of the grav- 
itational force and the external force F™: F™ = 
~f7™^(x)e^^(a;)i;^'y*=^ + F™. The backward Kolmogorov 
equation for the stochastic process described by Equ. (25) 
is defined by 



d_ 

a7 



ip{T, X, V) = Ap(^(T, X, V). 



(28) 
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The coresponding Fokker-Planck equation in phase space 
(general relativistic Kramers equation) within the frame 
of general relativity is defined by the adjoint of the op- 
erator Ap. Since the Laplace-Beltrami operator is self- 
adjoint = A+ it takes the form: 



97 



v'^'^div^{eM{x)<i>) - div^{F<^) + ^A,„$, (29) 



where A^ is the Laplace-Beltrami Operator in the hyper- 
bolic velocity space 



A„ = C?'"«_-^!_ + G""7^„T^ (30) 
1 (9 8 

The divergence operator in the position space is given by 

1 d 

div^{eM{x)<i>) = ^^(V5eAf(a;)*), (31) 
and in the velocity space by 



Gdv' 



(32) 



with G = det{Gij}, g — det{gij}. 

Equ.(29) is the diffusion equation in the parametriza- 
tion of the phase-space proper time within the frame 
of general relativity for the probability density func- 
tion $ — (/)(t;x, v) or the the transition probability 
<I>(x,v,T I Xo,Vo,0). 

For the solution of the relativistic diffusion equation 
it is convenient to introduce the hyperbolic coordinate 
system for the 4- velocity defined by = shasini9cosiy9, 

= shasini9sin(/3, — shacosi? and w° = cha. We 
denote the velocities in the non-Cartesian coordinates 
by — a,v'^ = 6,v'^ — (p ,a — 1,2,3. The metric in 
this coordinates are simply to calculate and are given by 
Gil — 1, G22 — sh'^a, G33 — sh'^asin^ and Gij = for 
i ^ j,G = sh^a sin^ d. With the given metric the Laplace 
Beltrami Operator Au takes the form 



52 „ , a 

A = — ^ + 2ctha— 



(33) 



1 



92 



(sha)2 \di} 



ctgt9— -I- — 



1 92 



(sin ^?) 2 5(^5 



and 



divjF$) ^ (sha)-2 ((sha)2F"$)- 
oa 

-(sha)^i(sin79)-i^ (sini9F''$) 

-(sha)-i(sin79)-i^(F'^$), 
dip 



(34) 



is the divergence operator. Here the force components in 
the hyperbolic coordinate system F" , F"^ , F"^ are related 
with F' by = (cha)~^ [sin 7? (cos (^F^ -I- sinipF'^) + 
cos'dF^], F^ = (sha)-i [cos ■& {cos ipF'^ + sin ipF'^) - 
sin'&F^] and F^ = (sha)-i (sin??)-i[- siniyj^i _^ 
cos ipF'^]. 

The relativistic diffusion equation (29) is parameter- 
ized in terms of the phase-space proper time t. But it is 
more convenient to parameterize the stochastic process 
alternatively in terms of the coordinate time because the 
gravitational fields and the external force fields are given 
in this parametrization. In general relativity the observer 
time with the infinitesimal element dx^ = e\f{x)v^^ dr is 
a function of the proper time r and the space and velocity 
variables. This general definition introduces difficulties in 
the diffusion formalism. In order to avoid such problem 
and to simplify the derivation we describe diffusion in a 
frame system in which the time-like components of the 
frames em{x) vanish: e'^{x) = 0. This condition can be 
achieved in general if we impose three frame subsidiary 
conditions and remove a part of the six frame arbitraries. 
In particular, the condition e^(x) = is introduced in 
the 1-1-3 spacetime slicing in the Arnovitt-Deser-Misner 
(ADM) formalism [35*] of the hamiltonian formulation of 
gravity. In the ADM treatment the space-time manifold 
is split into a one-parameter family of space-like hyper- 
surfaces I](i) parameterized by a time-like function t or 
as a foliation of the hypersurfaces t = const. 

With vanishing frame components e'^{x) = the 
proper time is given by dx^ — (5oo)'w'^dT where u° is 
defined by the relation (17). The diffusion equation in 
the parametrization of the observer time can be derived 
from the stochastic differential equation (25) using the 
mathematical theorem of random time change in stochas- 
tic differential equations (see e.g. [l^l, (Si)- The proper 
time T is related with x'^ by the random transformation 



N{s){v"{s))-^ds, 



(35) 



with ~ (goo)~2 , Since r depends only on the stochastic 
events earlier than x'^ this random time change is an 
adapted process and therefore the time change of an Ito 
integral is again an_Ito integral, but driven by a different 
Wiener process dW{x°) = dW{T)N-^v°)^/^ [s^. 
This rule for a random time change is valid within the Ito 
calculus. Using the transformation rules of an Ito integral 
into a Stratonovich integral and dr = {rf')~^Ndx^ the 
relativistic Langevin equation (25) can be rewritten in 



6 



the parametrization with a;" as follows 

dx'ix°) = e'M{x)v^'\v°)-^Ndx^, 



(36) 



D 



On-1/2 



{{v^)~^/'^)dx^ + F"(^>°)-liVdx^ 



d 

~^Zi^)Elodv-{x' 



Note that in Equ. (36 ) an additional drift term pro- 
portional to the diffusion constant D arise which origi- 
nates from the transformation rule from the Ito to the 
Stratonovich calculus for random time changes. The dif- 
fusion operator Kp^M) in the fiber bundle F{Ml) with 
coordinates r ={x*, w™, i?™} for the stochastic process 
described by Equ. (25) is defined as in Equ. (5) with hor- 
izontal vector fields Ha and Hq now given by 



E^{v°y^^^Ni 



d 



(37) 



-jZ{v)ElEl\v°)-^/'Ni 



dE" 



{-^S^'>ETEl^N{v")-'/^^{{v" 



rbc ™ pfc_^// ,0wl/2\l 



with the gravitational and the external force F™ = 
~n'i'^j^{x)e'^j{x)v'^v'^'' + F^. The operator Af(m) can be 
projected to the phase space with coordinates r ={x*,?;'} 
by Aj-(Af)/(r) = Ap/(x, v), where the diffusion genera- 
tor in the phase space Ap is given by 



(38) 



+elfix)v^'iv°)-'N^ + F'"(t;°)-i7V-^ 

The general relativistic Kramers equation in the 
parametrization of the observer time xq can be derived 
analogous as above and is written as follows 

= -v^'div,{eM{x)<^>) - c^^^;,(F$) + yA^*. 

(39) 

As seen in the limit of special relativity the left side 
and the first term o.r.s. of Equ. (39) can be identified with 
the covariant expression w'^S/c^a::^, while the other terms 
are identical with corresponding terms in Equ. (29). 

The equation (39) is the main result of the present 
paper and represents the generalization of the Kramers 



equation within the frame of general relativity in the 
parametrization of the observer time for the probability 
density function $ = (/)(xo;x, v) with the initial condi- 
tion (j>{xQ = 0;x, v) = 0q(x, v). The transition prob- 
ability is determined by the same equation but is de- 
fined by the initial condition $(x, v, | xo,vo,0) = 
{Gg)-^/^S{v^ - vl) d{v^ - vi)S{v^ - v^)S{xi ~ x'i)S{x2 - 
X2)S{x3 — X3). For an external electromagnetic field 
F^i, the normalized force in the moving frame is 
F™ = e0'^g^PFp,e''^{x)v^. 



4. COORDINATE TRANSFORMATIONS 

In the general relativistic framework, the invariance 
of the physical laws with respect of general coordinate 
transformations is one of the most fundamental property. 
In the frame of special relativity the probability density 
function (/)(r; x, v) in the phase space is Lorentz invariant; 
i.e. it fulfills the condition 



0'(r' 



(40) 



where the variables x' , v' are related with x, v by a 
Lorentz transformation. Note that in contrast the parti- 
cle density and the current density (i.e. the integrals over 
the velocities) transform like a four-vector. As shown in 
3l| the special relativistic equation (9) suffices this con- 
dition and is invariant with respect to a Lorentz trans- 
formation. In the general theory of relativity this invari- 
ance should be valid for general coordinate transforma- 
tions. However in the derivation of Equ. (39) using the 
parametrization of the observer time we have taken ad- 
vantage of the freedom in the choice of the orthonormal 
frame components and used the condition e^(a;) = or 
correspondingly a hypersurface foliation. In this frame a 
coordinate independent notion of time is demanded, and 
therefore we have to consider a foliation-preserving dif- 
feomorphism described by the general space coordinate 
transformations 



x"=f{x\t), 



(41) 



which preserve the hypersurface geometry. The hyper- 
surface foliation is also preserved under an arbitrary time 
rescaling 



t' = fit) 



(42) 



which do not affect the hypersurfaces. General covari- 
ance will then become a hidden feature similar as in the 
hamiltonian formulation of gravity in the ADM formal- 
ism, but the underlying invariance of general relativity 
is intact and general coordinate transformation still map 
solutions into solutions. 
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Let us study the transformation property of the diffu- 
sion equation (39) in the parametrization of the observer 
time. Since the moving 1-forms 9^{x) transform Uke a 
co-vector we find from the relation Equ. (16) that the 
vector components w*^ in the moving frame are invari- 
ant with respect to general coordinate transformations 



F . Therefore the operators divy and 
in the last two terms o.r.s. of Equ. (39) are also in- 
variant. From the transformation described by Equ. (41) 
we find with = x"^ 



dp d 
dx^ dx'i 



(43) 



The first term o.r.s. of Equ. (39) including the three- 
dimensional divx operator with respect to the 3-geometry 
of the hypersurface is intrinsically defined by the hyper- 
surface and therefore invariant under the hypersurface- 
preserving transformation (41). This can be simply 
proven by the transformation property of the intrinsic co- 
variant differentiation '^^^ AiF^ on the hypersurface given 
by (3)A'jF" = {dx" /dx''){dx^ /dx'3)^^^ AiF'' . For the 
divergence operator diVx^Fcf)) — '^■'Aj(F-'0) this yields 
with [dx'' / dx^){dx^ / dx'') = 5^ the relation (3)A'jF'J 
=(^)AjF-' or diVx'{F'<p) — diVxiFcf)). The time derivative 
on the left-hand side transforms under Equ. (41) as 



dx 



Ik 



(44) 



ok ^ df" dx^ 



The transformed relativistic Kramers equation (39) 
takes therefore the form : 



5. DIFFUSION IN THE EXPANDING UNIVERSE 

Let us discuss the above given relativistic diffusion 
equation (39) in gravitational fields for the example of 
the expanding universe. The metric that characterize the 
expanding spatial homogenous and isotropic universe can 
be described by the Robertson- Walker metric which has 
the form 

ds' = [cAtf - P^[{dx^f + {dx^f + {dx^fl (47) 

where e may assume the values 0, lor -1 for a flat, closed 
or open universe, respectively, r^ = {x^Y + {x^Y + {^^Y 
and R{t) is the cosmic scale factor. In the following we 
restrict ourself to the spatial flat case with e = 0. With 
the metric coefficients goo — ^^,9Qi = 0, Qa = —R{t)'^Sij 
we find for the Christoffel coefficients [l^, 



TjTc = 0> r°fe = R-^djkXok = i^R'^^ik- (48) 



dR 



The moving frames related with the metric in Equ. (47) 
are given by 

e'; = i?(x°)-i<5°,e° = l^ei = R{x')5iX = 1. (49) 

Using Equ. (48) and (49) we find for the spin connection 
coefficients 



O"^ — no" — A° — n 



(50) 



The gravitational force F^ — —Q'^^^g{x)e^{x)u^u^ is 
given by 



= -«'^fc(e:,(x)$) (45) 
-diw„,(F$) + yA„,$, 

where 

div^,{e'^{x')<^) = -L_^(^4(x')$), (46) 

and g' = gj^ where J is the Jacobian matrix of the trans- 
formation: J — det{^^}. A time rescaling Equ. (42) do 
not change the left-hand side of Equ. (39) because the en- 
tity N is transformed like N' — N{dt/dt') and therefore 
iN)-^d/dt = {N')-^d/dt'. 

Note that transformation properties as found in Equ. 
(45) with the replacement of the time derivative by the 
left-side in Equ. (45) is a general feature of the evolu- 
tion formalism in the 3-1-1 spacetime decomposition in 
the general relativity theory (see e.g. [s^l). 



- -d^^^'"^"^"- (51) 

Let us discuss the diffusion of massive particles in 
the expanding universe. Introducing the hyperbolic co- 
ordinate system for the 4-velocity, defined by = 
shasiniJcosiyj, v? = shasini?sin(/j, = shacos?? and 
vP = cha we find for the gravitational force V^{ti) — 
~H{rj)sha, F^ ^ F^ = where H{r]) = |f i?"^ is the 
Hubble constant and r] = xq = rch. The random im- 
pact of surrounding particles generally cause two kind of 
effects: they act as a random driving force leading to a 
random motion and they give rise to a frictional force. In 
the non-relativistic theory the friction force is given by 
fp — —vmv^ where v is the friction coefficient and are 
the components of the non-relativistic velocity. The rel- 
ativistic generalization of the friction force requires the 
introduction of a friction tensor i/^ similar to the pressure 
tensor in the relativity theory 01; 0. The friction force 
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is expressed as Fp = v'^J\u°' — U"], where C/" is the 4- 
velocity of the heat bath. For an isotropic homogeneous 
heat bath the friction tensor is given by 



(52) 



with V denoting the scalar friction coefficient measured in 
the rest frame of the particles. In the laboratory frame 
the heat bath is at rest described by t/° = (1,0,0,0). 
Therefore the friction force is given by Fp = —vu''u^^ or 
in hyperbolic coordinates 



-i/sha, Fp 



^ F 



0. 



(53) 



In the expanding universe the coefficients i^(?7) and 
D{'q) depend on time and it is convenient to use the dif- 
fusion equation (39) in term of the parametrization with 
the observer time 77. In the case of a spatial homogenous 
and isotropic solution in Equ.(39) the spatial derivatives 
vanish. Substituting the ansatz cj) ~ (l)'l[(a,'d,ip,r]) — 
gj{a,r^)Yl,{e,ip) with Yi,{e,^) =P5[,(z?)e'^'^ as the 
spherical harmonics and the associated Legendre func- 
tions into Equ. (39) the following equation can 
be derived for gj = gj{a,r]) : 



cha-^gj 
or] 



^ 2 ^da' 



2cthQ! 



d J{J + 1) , 



dc 



sh^ a 



-x(?7)(shQ:) ^— sha)^]gj. 



(54) 



with xiv) = §f + Here the discrete index J 

takes the values J = 0,1,2,... and M = —J,—J + 
1, ... 0, 1, ... J. For the fundamental solution J = we 
find the diffusion equation 

cha-go = [^(a^ + 2ctha-)-f 

-|-x(?7)(sha)-2Ash3a]5o. (55) 
Let us discuss the solution of this equation within a cer- 
tain range of validity substituting the special ansatz 

4>{a, r/) — C exp{/3(ry) — 7(77) cha} (56) 
into Equ. (55) which yields the relation 

(57) 



^ -^cha = x(?7)[3 — 7s/iq</iq] 
or] or] 



+ —^(-37 + -f^shatha). 

In general, there do not exist functions I3{r]) and 7(77) 
which solve this equation, but if we restrict ourself to the 



ultrarelativistic case a ^ 1 we find as equations for the 
coefficients /3(ry) and 7(7/): 



(58) 



The constant f3 is included into the normalization 
of (j){a,r]). An analytical solution of Equ. (58) for 7(77) 
can be found by a transformation of the variable 7(77) = 
{Y(rj))~^ which yield the solution 



^°i?(77)' +2i?(7y)' 



Hv) I p''(*')i?(t')i:»(t')di' 



(59) 



with h{t) = / v(t')dt'. 

to 

By comparing the Jiittner distribution with the solu- 
tion ansatz (56) we can introduce a time-depending tem- 
perature of the expanding universe 



7(77) — m(?{kT{r])) ^ . 



(60) 



Since the diffusion and drift constants D and v are de- 
termined by the scattering processes of the particles in 
the system described by different physical parameters; 
in particular by the temperature, the solution (59) has 
for temperature depending friction and diffusion coeffi- 
cients the meaning of an integral equation. However, we 
have to taken into account that the diffusion and friction 
coefficients are related each others by the fluctuation- 
dissipation theorem. This relation is well known in the 
non-relativistic case where the viscous friction coefficient 
v oi a Brownian particle must be related to the diffusion 
constant D of the particles by the Einstein relation 



D = 



2vkT 



(61) 



The nature of the random force is independent on the 
presence of the gravitational field. Therefore in the rela- 
tivistic case the stationary solution of the Equ. (55) for 
^Rr^ — ^ must coincide with the Jiittner distribution. 
From the recently derived special relativistic diffusion 
equation 's^] it was shown that the Jiittner distribution 
arise as the stationary solution of the special relativis- 
tic diffusion equation (9) if the Einstein relation is not 
only valid in the non-relativistic case, but also in the rel- 
ativistic regime. Substituting the relation (61) with (60) 
into Equ. (58) one can see, that the effect of diffusion is 
canceled by the viscosity. Then from Equ. (58) we obtain 



liv) = 7o 



Ro 



(62) 
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This relation is identical with the result in the kinetic the- 
ory for the expanding universe [s^ l , (35| . In the radiation 
dominated period in a flat cosmos we have R{i]) ^ y/rj 
and in the matter dominated period R{r]) ^ rf^^ and 
therefore we find for the ultra-relativistic case T ~ 77"^/^ 
or T r\~'^l^ ^ respectively for the corresponding periods. 

Note that the solution (56) do not satisfy physically 
determined initial conditions; this solution describes the 
asymptotic quasi-static regime and is valid only after a 
certain time when the system is already in the equilib- 
rium. As shown above, in this case the diffusion is com- 
pensated by the drift process, and the result that follows 
is consistent with the equilibrium state in kinetic theory. 
On the other hand in the opposite transient case up to a 
certain time 77 after the initial time r^g one can neglect in 
Equ.(55) the last term proportional to x('7)- Then in the 
ultra-relativistic case the fundamental solution J = is 
determined by the equation 



5^ 1.. 5^ o^N^ 



(63) 



where the new variables ^ — J D{t)dt and ^ = exp(a) are 



introduced. Using the Laplace transformation <1'(<;,^) 

00 

/ ^(A,^) exp(— A<;)(i<; we find the solution 




$(A,0=rVi(2v/|A), 



(64) 



where Ji is the first order Bessel function. The eigen- 
functions $(A,^) satisfy the relations of orthogonality. 
Therefore the transition probability is determined by 



<f(e,? I Co 0) = j $(A,^ )$*(A,eo)exp(-AOdA. (65) 


Substituting Equ.(64) into Equ. (65) we find 



-)exp(- 



(66) 



where C is the normalization constant and Ii{x) — 
—iJi(ix). In the general case of arbitrary time the so- 
lution of Equ. (54) can be obtained by numerical meth- 
ods. But for the study of this problem under the con- 
ditions of the earliest epoch of the universe we need a 
realistic microscopic model for the viscosity in the non- 
equilibrium epoch in a plasma of relativistic particles, 
including quarks, leptons, gauge and Higgs bosons. A 
detailed discussion of this issue is beyond the scope of 
the present paper. Corresponding cosmological estima- 
tions the universe may not have been in thermal equi- 
librium during its earliest epoch in a time range earlier 



than about 10~^®s after the big bang or temperatures 
greater than IQ^^GeV [s^. Standard phenomenological 
inflationary cosmological models relate this epoch with 
symmetry breaking phase transitions. 

Let us finally briefly discuss the non-relativistic limit 
a <C l.From Equ. (57) we find under this condition: 



97 
dr] 

dr] 



(67) 



If we use the Einstein relation (61) the diffusion term 
is again canceled by the friction and from Equ. (67) the 
solution 7(77) = 7g i?^(77)/i?g follows. The same solution 
is obtained in the kinetic theory for a non-relativistic 
gas. In the non-relativistic case the temperature of the 
equilibrium distribution therefore depends on the cosmic 
scale factor R{r]) like T ^ R^^{rj). Thus, both in the rela- 
tivistic and in the non-relativistic case the validity of the 
general fluctuation-dissipation theorem with the Einstein 
relation (61) ensures that the result of the kinetic theory 
derived from the vanishing of the collision integral in the 
Boltzmann equation is consistent with the here derived 
probabilistic general relativistic diffusion theory in the 
quasi-steady state regime or in thermodynamic equilib- 
rium. 



6. CONCLUSIONS 

In conclusion, a theory of Markovian diffusion pro- 
cesses within the framework of the general theory of 
relativity is formulated. In the derivation of the ba- 
sic relativistic diffusion equation the mathematical cal- 
culus of stochastic differential equations on Riemannian 
manifolds is used, which here is modified for the descrip- 
tion of diffusion in the phase space of Pseudo-Riemanian 
manifolds with an indefinite metric by using orthonor- 
mal frame vectors both in the position and in the velocity 
space. A generalized Langevin equation in the fiber space 
of position, velocity and orthonormal velocity frames is 
defined and the generalized Kramers equation within the 
framework of general relativity is derived both in the 
parametrization of the phase-space proper time and the 
observer time. The transformation of the obtained diffu- 
sion equation under hypersurface-preserving coordinate 
transformations is studied and diffusion in the expand- 
ing universe is discussed. It is shown that the validity 
of the fluctuation-dissipation theorem in the relativistic 
case ensures that in the quasi-steady state regime the re- 
sult of the derived diffusion equation is consistent with 
the kinetic theory in thermodynamic equilibrium. Be- 
sides a transient analytical solution valid for small times 
has been derived. 
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